SPECIAL VALUES OF PARTIAL ZETA FUNCTIONS OF REAL 
QUADRATIC FIELDS AT NONPOSLITVE INTEGERS AND 
EULER-MACLAURIN FORMULA 

BYUNGHEUP JUN AND JUNGYUN LEE 

q_) ■ Abstract. We compute the special values at nonpositive integers of 

C/3 . the partial zeta function of an ideal of a real quadratic field applying 

an asymptotic version of Euler-Maclaurin formula to the lattice cone 
associated to the ideal considered. The Euler-Maclaurin formula in- 
volved is obtained by applying the Todd series of differential opera- 
tors to an integral of a small perturbation of the cone. The additive 
property of Todd series w.r.t. the cone decomposition enables us to 
express the partial zeta values in terms of the continued fraction of 
the reduced element of the ideal. The expression obtained uses the 
positive continued fraction which yields a virtual decomposition of 
the cone. 

We apply the expression to some indexed families of real quadratic 
fields satisfying certain condition on the shape of the continued frac- 
tions. The families considered include those appeared in II20II and 
II21II as well as the Richaud-Degert types. We show that the partial 
zeta values at a given nonpositive integer — k in the family indexed 
by n is a polynomial of n. 

Finally, we compute explicitly the polynomials producing the par- 
tial zeta values at s = — fc for small k of some chosen families and 
compare these with some previously known results. 
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1. Introduction 

Let K be a number field of the extension degree [K : Q] = r 1 + 2r 2 , 
where r 1 and r 2 denote respectively the number of real and complex 
embeddings of K. The Dedekind zeta function 

T r 1 

11 1 -Afp- S 

p :prime ideal in K r 

is encoded with many interesting arithmetic properties of K. In particu- 
lar, the residue at s = 1 is associated to the class number h K of K by the 
class number formula: 

2^(2n)^R K h K 

Res s=1 £ K (s) = ^— , 

^kV\D k \ 

where R K is the regulator, co K is the number of roots of 1 in if and D K 
is the discriminant. This has been the starting point of most studies of 
class numbers. 

The simplest is the case of imaginary quadratic fields where the reg- 
ulator appears to be trivial. In 111311 . Gauss listed 9 imaginary quadratic 
fields of class number 1 and conjectured that the list is complete. Later 
on this had been studied through 20th century and is now quite well 
understood and solved by works of Heegner, Stark, Goldfeld and several 
others(eg. [ffl], EzD, D3TD, ESD, CUD, ESD, EH and 11321). 

The case of real quadratic fields is more complicated due to the pres- 
ence of nontrivial regulator. It is also conjectured by Gauss that there 
are infinitely many real quadratic fields of class number one. But since 
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the regulator is far from being controlled in relation to the discriminant, 
and there has been no essential progress to the proof of the conjecture. 

Instead of treating the whole real quadratic fields, people considered 
some families of real quadratic fields where the regulators are controlled 
in relation to the discriminant. The most well-known family of this kind 
is the Richaud-Degert type: A Richaud-Degert type is defined by 

d(n) = n 2 ±r 

for r\4n and n < r ^ n. For r fixed as above, the family {K n } of real 
quadratic fields is called R-D type. In this case, we have a bound of the 
regulator R Kn : 

R Kn <3\og^D^ n 

As in imaginary quadratic case, a well-known estimation of Siegel L(l, % D ) ~ 
\D\~ £ together with the class number formula implies that there are only 
finitely many R-D type fields of class number one. Assuming the gener- 
alized Riemann hypothesis, the class number one problems have been 
solved for many subfamilies in R-D type. 

It is quite recent that Biro first obtained an Riemann hypothesis free 
answer to the class number one problem for the families K n = Q(V n 2 + 4) 
and K n = Q(\/4n 2 + 1) in a series of papers CQ2D , O)- He investigated 
the behavior of the special values of the partial Hecke L-functions at 
5 = in the family. The partial Hecke L-function of an ideal a is defined 
for a ray class character % as 

He discovered that the special values behave in a packet of linear forms 
whose coefficients are easily computed for the family (K n , Kn ,x„ '■= X ° 
N k j<q) for a Dirichlet character %■ This property is named the linearity. 

Inspired by Biro's pioneering work, in H5J, fl6|], II25II and II26II the 
linearity is observed for more general families of Richaud-Degert types 
and the class number one and two problems have been answered for 
these. 

In | |20 || , we found a sufficient condition to yield the linearity of the 
Hecke L- values at s = 0. Namely, for families of integral ideals {b n in 
K n } such that b^ 1 = [1, o>(n)] := 1A + Z(o(n), where co(n) has purely 
periodic positive continued fraction expansion of a fixed period r 

w(n) = [[ao(n)» ai («).••• .flr-i(n)]] 

and N(b n )N(xco(n) + y) = b (n)x 2 + b 1 (n)xy + b 2 (n)y 2 for a^n) and 
b;(n) being integer coefficient linear forms in n. In this setting we have, 
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for n = qk+r with < r < q and for a Dirichlet character % of conductor 
q, the L -value at s = 

L Kn (0, X oN Kn ,b n ) = J-j (A z (r)fc+fl,(r)) 

with A„(r), B x (r) e Z[x] where Z[x] denotes the extension of Z by the 
values of % . 

In | |21 || we obtained a higher degree generalization of the linearity for 
ray class partial zeta values. Let co(n) be the Gauss' reduced element of 
b n . If we allow the coefficient a ; (n) of the continued fraction of co(n) to 
be polynomial of degree d, then the partial zeta value at s = of a mod- 
q ray class ideal (C + Dco(n))b n in the class of b n is a quasi-polynomial 
in n: 

C q (0, (C + Dco(n))b„) = ^ (A (r) +A a (r)k + • • • + A d (r)k d ) 

with A; (r) e Z(for precise definition, we refer the reader to loc.cit). In 
particular, if we take d = 1 and sum the ray class zeta values twisted 
by x n , one can recover the linearity of the partial Hecke L- values. For 
d > 1, the same process concludes the polynomial behavior of the partial 
Hecke values at s = 0. 

The purpose of this article is to generalize our earlier work to spe- 
cial values at every nonpositive integer of the ideal class partial zeta 
functions under the same assumption for the family (K n , b n ). We assume 
again b" 1 = [1, co(n)] where co(n) has purely periodic continued fraction 
expansion of fixed period r: 

a>(n) = [[ a oO)> a iO)>-" , a r-i(n)]] 

such that and N(b n )N(xco(n) + y) = b (n)x 2 + b x {n)xy + b 2 (n)y 2 for 
a f (n) and b^n) being integer coefficient polynomials. For the next two 
theorems, let I be the even period of co(n) (hence independent of n and 
£ = 2r(reps. r) if r is odd(reps. even)). 
Our main result in this paper is as follows: 

Theorem 1.1. Let N be a fixed subset ofN. Suppose (K n ,b n ) satisfies the 
above condition for every n £ N. Then the special value of the partial zeta 
function of b n at s = —kfor k = 0, 1, 2, • • is given by a polynomial in n: 



Ck„(-M„) =A +A l n+A 2 n 2 + ...+A m n' 
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of degree bounded by m = kC + D with the coefficients A f e -rZ, where C, 
D and C fc are given as /oHows: 

C = 2 deg ct£ _! + deg b x 



and 

C k = LCM of 



D = max {deg a,} 

OsgisSr-1 



^i + Afc+l-i , ^2fc+2 f 2fc N 



{ i + 1 2k+l—i 
the denominators of — and 



(i + l)(2fc + l-i) (2fc + 2)(2fc + l) V 1 

Our main theorem is a direct consequence of the following estimation 
of the partial zeta values of an ideal b in a real quadratic field K. 

Let b be an integral ideal such that b _1 = [1, co] for co > 1 and < 
co' < 1, Let a;,/?; are coordinates of some lattice vectors determined by 
(the continued fraction of) co. In particular, a^co + /3^_ 1 1 is the totally 
positive fundamental unit of K = K n . See Secj6]for details. As usual, B t 
denotes the i-th Bernoulli number. 

Theorem 1.2. Let b be an ideal of a real quadratic field K such that b _1 = 
[1, co] where co = [[a , a 1 ,..., a r _ l ]]. Then we have 

a-Kb) = 

i-i 



i=0 



R 

■02^+2 



+ ( 2 fc"+2)! S^ -1 )^-^^!* d h 2 )Q( a i-2hi + a t h 2 , p t _ 2 h x + fah 2 ) k 

In the above, L fc and R k are the homogeneous polynomials of degree 
2k: 

2fc+l R R 

(1.1) L k (X,Y) = V _i x f-i y2 fc-£ + i 

fcV 7 £J i! (2fc + 2 — £)! 

(1.2) R*(X,Y) = X 2fc +X 2fc_1 Y+-+Y 2fc . 

It is not surprising that this behavior of the partial zeta or L-values is 
related to the pattern of the continued fractions in the family if we note 
that the Shintani cone decomposition arises in relation to the continued 
fraction. The significance of this fact lies on that for real quadratic fields, 
the regulator is controlled not only by the discriminant but also by the 
period r of the positive continued fraction of the reduced element co in 
K. This is due to the following well-known upper bound: 

R K ^r\og^D~ K . 
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Study of special values of zeta or L-functions goes back to Euler. Euler 
evaluated £(— fc) for k = 0,1,2,..., by using Euler-MacLaurin summa- 
tion formula(cf. []8l]). Later on, Siegel computed the values at non- 
positive integers of £ K (b,f,s) the ray class partial zeta function for an 
ideal b in a totally real number field K w.r.t. a conductor f based on 
the theory of modular forms rf|30H"). Shintani established a combinato- 
rial description of the zeta values at nonpositive integers ([29]). Beside 
the complicated contour integral, Shintani's method is a reminiscence 
of Euler's. Similar approach was taken independently by Zagier in his 
evaluation of the partial zeta functions of real quadratic fields at non- 
positive integers (Q34D). Actually, the Shintani's method has a strength 
over Siegel's that it is ready to use in p-adic interpolation in case of to- 
tally real fields via Cartier duality. This view was clarified by Katz in 
Q24fl. 

Our evaluation is along with the line of Shintani and Zagier. We will 
apply a version of Euler-Maclaurin summation formula due to Karshon- 
Sternberg-Weitsmanf ||22|| 1 . In loc.cit, they made a version of Euler- 
Maclaurin formula taking care of the remainder term, so that one can 
apply this to expand asymptotically a function given as summation of ex- 
ponentials. Since one side of the Euler-Maclaurin formula is application 
of appropriate version of Todd differential operator, the decomposition 
of Shintani cone is reflected additively due to the additivity of the Todd 
series under cone decomposition(See SS J4.4D . 

Similar computation was done by Garoufalidis-Pommersheimr||12||'). 
They applied the Euler-Maclaurin summation formula of Brion-Vergne ( [ffl ) 
to obtain the asymptotic expansion. Brion-Vergne's formula is exact sum- 
mation on the lattice points inside a simple polytope valid for polynomi- 
als or polynomials in exponentials of linear forms. They took the Shin- 
tani cone as the cone over a lattice polytope and varied the size inside 
the cone. In their treatment, the exact and the error terms are consid- 
ered separately. The formula of Brion-Vergne is applied to the exact term 
and the error term is shown to be appropriately bounded. 

Our method differs from that of Garoufalidis and Pommersheim in 
two directions: First, we used positive continued fraction while they 
took negative continued fractions. Basically, via the transition formula 
between positive and negative continued fractions, they contain more or 
less the same information of the ideal. But as is pointed at the beginning, 
it is important to note that the period of the two continued fractions have 
no control on each other. As the regulator is concerned, it is better to 
express the zeta values using the terms of the positive continued frac- 
tion. Nevertheless, our earlier work has been made via translation of 
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the terms of the positive continued fraction into those of negative con- 
tinued fraction. Hence the direct use of the positive continued fraction 
significantly reduces the amount of the computations needed. While the 
negative continued fraction yields an actual cone decomposition, the 
positive continued fraction gives rise to a virtual cone decomposition. In 
fact, the cone decomposition appeared at first to be a fan in toric geom- 
etry where no virtual decomposition is allowed to define a toric variety. 
The Todd additivity can be simply extended to virtual decompositions 
if we take care of the orientation of the cone and take it as the sign of 
the Todd series. Second, we make a direct use of the Euler-Maclaurin 
formula of Karshon-Sternberg-Weitsman. Since the exponentials in our 
case is of Schwarz class toward the infinity of the cone to evaluate. Thus 
both sides of the Euler-Maclaurin formula make sense and we have more 
transparent proof. In addition, Karshon-Sternberg-Weitman's version of 
Euler-Maclaurin has advantage over Brion-Vergne's in that no limitation 
on the function to integrate while Brion-Vergne's, since the former can 
be applied to wider range of functions. 

Again, some part of the computation is similar to what had been done 
by ZagierC [13410 . He obtained the partial zeta values at non-positive 
integers again by decomposition of the underlying cone of the zeta sum- 
mation according to the negative continued fraction together with the 
Euler-Maclaurin formula. One should note that this decomposition ap- 
pears in the dual side while the additive decomposition of the Todd dif- 
ferential operator is taken in this paper and H12 || . It was already pointed 
out in | |12 || as "M-additivity" to "N-additivity". One could ask direct 
application of Zagier's method with the cone decomposition from the 
positive continued fraction. Unfortunately, in this setting, the both sides 
of Euler-Maclaurin formula don't make sense but need certain renormal- 
ization process to avoid operations with infinity. In our setting, this is no 
problem as the domain of integration remains the same while the cone 
decomposition is reflected to the Todd differential operator. 

The plan of this paper is as follows: First, we rewrite the partial zeta 
function of an ideal as a zeta function of a quadratic form weighted by a 
fundamental lattice cone of the Shintani decomposition(Sec. [2]). We re- 
call a standard asymptotic method to evaluate the zeta values at nonpos- 
itive integers and rebuild a version of Euler-Maclaurin formula (Sec|3}- 
[4]) . Then we apply this Euler-Maclaurin formula to obtain an expression 
of the zeta values and another expression after the cone decomposition 
arising from the positive continued fractions (Sec |5][7]). The partial zeta 
values at 5 = 0,-1,-2 are explicitly computed using our method and 
compared with previously known results (Sec JU). Secj9j which is tech- 
nical and similar to the computation by Zagier( [j34l "), is devoted to the 
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proof of vanishing of a part skipped in the previous sections. Finally, 
we apply this to some families of real quadratic fields to prove our main 
theorem ("Thm J 1 . ll) and the polynomials are explicitly computed out for 
some families and for some small kfSec JTOl) . 

2. Partial zeta function of real quadratic fields 

2.1. Partial zeta function. Let K be a real quadratic field and b be an 
ideal. Through out this article, by partial zeta function, we mean the 
partial zeta function of an ideal class in narrow sense. The partial zeta 
function of an ideal b is defined as 

f(s,&):= J] N(ay 

adntegral 

where a ~ b means b = aa for totally positive a in K. This infinite series 
defines a holomorphic function in the region Re(s) > 1 of the complex 
plane and has a meromophic continuation to the entire complex plane. 
Since for an integral ideal a in the narrow class of b there exists totally 
positive element a e b^ 1 such that a = ab and vice versa, we can write 
again 

C(s,b)= J] N(ab)- S , 

[a]e(b- 1 ) + /E+ 

where (b _1 ) + denotes the set of totally positive elements of b _1 and 
E + = denotes the group of totally positive units of K. Now we are 
going to describe the summation as taken inside the Minkowski space of 
K. Let (t l5 1 2 ) be two real embeddings of K. Let us denote the Minkowski 
space of K by 

K R =K® q R = K ti xK l2 

Then one can identify an ideal c with a lattice of K R given by its image 
under the diagonal embedding of K into K R : 

l = (i l3 t 2 ) :K^K M , 0(a) ^ (^(a), t 2 (a)) 

This is a full lattice in the Minkowski space. 

E£ acts on the 1st quadrant of K R by coordinate-wise multiplication 
after the diagonal embedding. Let e be the totally positive fundamental 
unit of K. A fundamental domain of this action is given as a half-open 
cone F K of K R with basis {t(l), i(e)}: 

F K = {xl(1) +yt(e) e JR 2 |x ^ 0,y > 0}. 

For an ideal o of K or a lattice A of K m , we denote its intersection with 
F K by F K (a) or F K (A), respectively. 
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Figure 1. A fundamental cone of E, -action 



For [a] e (b 1 ) + /£ + , there is a unique representative a chosen in 
F K (b _1 ). Thus we have 

£(s,b) = J] N(ab)- S . 



2.2. Zeta function of 2-dimensional cones. Consider the standard lat- 
tice M := 1? in R 2 . Let Q(x,y) = ax 2 + bxy + cy 2 be a quadratic form. 
For two linearly independent vectors v 1 ,v 2 , let cr(v l3 v 2 ) be the cone in 
M 2 as the convex hull of the two rays M + v 1; IR + v 2 : 

°"0n v 2) : = + X 2 V 2 I x i > for i = 1,2}. 

For simplicity, we write a instead of cr(v 1 ,v 2 ) if v 1; v 2 is clear from the 
context. Following our convention on cones, the origin is not contained 
in a. 

Define a weight function wt a with respect to a as follows: 

1 lEint(cr) 
(2.1) wt<J (t) = \\ le 8(a)- (0,0) 

otherwise 

This strange weight is justified via identification of the partial zeta 
function with the zeta function of a lattice cone that will be defined 
soon below. The partial zeta function is a sum over the points of F K (b _1 ). 
Since the two edges of F K are related by the multiplication of the totally 
positive unit, the summands over both edges coincide. When we take F K 
as half-open cone, this repetition is automatically removed. Equivalently, 
we may apply this weight function so that the total contribution over an 
orbit equals 1. The choice of assigning 1/2 to each edge will be found 
useful when we apply the Euler-Maclaurin formula to the cone. 
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For a and a quadratic form Q(— ) satisfying Q(v) > for v e cr, we 
define a zeta function as the following series on Re(s) > 1: 



2.3. Comparison of zeta functions. One can choose b as an integral 
ideal in the same class such that b -1 = [1, to] for [l,to] being the free 
Z-module generated by 1 and to. Taking i(l),i(a>) as basis of K R , we 
have trivialization 



Here, we fix the order of the basis such that x + y to reads (y, x) in R 2 . 

From the reduction theory of quadratic forms, we have a privileged 
choice of to such that 



Then the totally positive fundamental unit e belongs to b 1 and e = 
p + qto for a pair (p,q) of relatively prime positive integers. 

Let cr be lattice cone generated by (0, 1) and (q,p), which corresponds 
to F K . One should be aware that this identification depends on b and the 
choice of to. Then we have for an integral ideal a = ab with a totally 
positive in K 



for a = mto + n e F K (b x ). Thus we have the following identification of 
zeta functions: 

Lemma 2.1. Let Q(m, n) = N(mto + n)N(b) and cr be a cone defined as 
above. Then we have 



3. Euler-Maclaurin formula and Zagier's asymptotics 

In the previous section, we have identified the partial zeta function 
of an ideal as a zeta function of a quadratic form defined over a lattice 
cone. 

To evaluate the values at nonpositive integers, we will apply Zagier's 
asymptotic method to the exponential series associated to the zeta func- 
tion of quadratic form running over the lattice points of the considered 
cone. The coefficients of the asymptotic expansion of the exponential 
series which will be obtained via Euler-Maclaurin formula are the zeta 
values at nonpositive integers up to some simple factors. 

We first recall the asymptotic method of Zagier then state the appro- 
priate Euler-Maclaurin formula for our case. 




K R ~ R 2 and ifV 1 ) M = I?. 



Ll(o>) > 1, — 1 < L 2 (to) < 0. 



AT (a) = N(mto + n)N(b) = Q(m, n), 



C(s,b) = ( Q (s,cr). 
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3.1. Zagier's asymptotic method. For a Dirichelet series of the follow- 
ing form 

Us):=2p W cR+ - A ^ go, 

X 

if meromophically continued to the entire complex plane, there is a fairly 
standard approach to evaluate the values at nonpositive integers. 

If Yjx a A e At nas asymptotic expansion Xl/l-i c ; f! at t = 0, then C A (s) 
has meromorphic continuation to entire complex plane and 

Ca(-h) = (-l) n n!c n 

for a nonnegative integer n(See Prop. 2 in [34]). 

Let Q(— ) be a quadratic form and cr be a lattice cone such that Q\ a is 
positive. The Dirichlet series defining the zeta function of (Q, cr) yields 
the following exponential sum 

(3.1) E(t,Q,cr):= J] wt a (£)e-^, 

iear\M 

where wt(-) is the weight function denned for cr in Sec. 2. 

The asymptotic expansion of the above exponential series will be com- 
puted after we state the appropriate version of Euler-Maclaurin formula 
with remainder for cr and the weight in consideration. 



4. Euler-Maclaurin formula for 2-d cones 

4.1. Twisted Todd and L-series. Let A be an iV-th root of 1. We define 
a A-twist of the classical Todd series. 

Todd^(S) = ^ 

When A = 1, this is the classical Todd series. This version of Todd series 
is used in [4] where they take the sum of the values of a function at the 
lattice points (strictly) inside of a simple lattice polytope. 

As we use weight 1/2 on the boundary rays of a cone, another variant 
of the Todd series with A-twist is defined as follows: for A a root of unity, 

Sl + Ae~ A S S 

L x S) = = 

V ; 2 1-Ae- A 1-Ae- S 2' 

We call this the A-twisted L-series. This fits well to the case when we put 
the weight (l/2) cod on a face, where 'cod' denotes the codimension of the 
face. This is used in | |22 || in their version of Euler-Maclaurin formula. 
For A = 1, L 1 (S) is nothing but the even part of Todd(S). Similarly to 
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the case of Todd series, the series expansion at S = of I A (S) is given 
as 

1 A 

L\S) = (-- — )S + Q 2 ,MS 2 + Q 3 ,a(0)S 3 + • • • + Q fc;A (0)S fc + • • • , 

where Q t A (x) is a (generalized) function of period N for i = 0, 1, 2, • • • 
defined as follows. 
For m = 0, 

For m > 1, Q m x(x) is an indefinite integral of Q m -ix( x ) with an integral 
constant fixed by the boundary value condition 

rN 

(4.2) Q mjA (x)dx = Q m+u (iV)-Q m+ljA (0) = 0. 

Jo 

Thus we have 

d [ N 
ax ' Jo 
Note that we are taking these Q A for m ^ as distributions on 
L*([0, +oo)). Q 1A is continuous and Q m A is C m_1 -function for m ^ 2. 
These generalize the periodic Bernoulli functions appearing in some lit- 
eratures on analytic continuation of the Riemann zeta function using the 
Euler-Maclaurin formula (eg. H8D ) . 

4.2. Todd series of 2 -dimensional cone. Let M = 1? a M 2 be a fixed 
lattice. Recall that a lattice cone is the convex hull of two rays generated 
by lattice vector. We may assume the generating vectors of a cone are 
primitive (i.e. not a multiple of other lattice vector in the same ray). For 
two linearly independent primitive lattice vectors v l3 v 2 , let <j{v x , v 2 ) be 
the cone generated by v a and v 2 . When v l3 v 2 are clear from the context, 
we will simply write cr intend of a(v 1 ,v 2 ). When there appear several 

cones, they will be denoted by a, t . . . or o~ l5 <r 2 , cr 3 , Since we will 

be concerned with surface integral over a 2-dimensional cone the order 
of the basis vectors (ie. the orientation of the cone) is important. So 
cj(v 1 ,v 2 ) is never equal to cr(v 2 ,v 1 ). Taking v l3 v 2 as column vectors in 
Z 2 , we associate a nonsingular (2 x 2)-matrix 

K = 0i, v 2 ) 

to a lattice cone cr = <r(v 1 ,v 2 ). Conversely, if a (2 x 2) -nonsingular 
matrix A with integer coefficient has column vectors v l3 v 2 which are 
primitive, we can associate a unique lattice cone. A cone is said to be 
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nonsingular if the matrix is in GI 2 (Z). Equivalently, cr is nonsingular iff 
det(A ff ) = +1. 

Remark 4.1. In literatures on polytopes or toric geometry, a cone is said 
to be simple if it is generated by n-linearly independent rays in M". In this 
article, as we are considering only 2 dimensional cones, every cone is simple 
unless degenerate. 

Let M G be the sublattice of M generated by v l ,v 2 and r ff be M/M a . 
An element g e r CT can be written as 

g = a aA (g)v 1 + a a2 {g)v 2 

for rational numbers a t7 l (g), a a2 (g) modulo Z. This is given ambigu- 
ously but yields two well-defined characters 

X*,i--g~e 2ataaAs) , fori = 1,2. 
The Todd power series for a cone cr is defined as 

Todd (J (x 1 ,x 2 ) := J] Todd^^^OTodd^^^). 

Similarly, we define the L-series for cr as 

L a {x 1 ,x 3 ):=%L*'.M( Xl ) L z.ato( x j m 

ger CT 

These are used in the Euler-Maclaurin formula of H and ||22||. re- 
spectively. For a 2-dim cone cr, the Todd and the L-series are related in 
the following manner: 

Lemma 4.2. 

|r ff | 

K{xi,x 2 ) = Todd <J (x 1 ,x 2 ) even - —x±x 2 , 

where Todd cr (x 1 ,x 2 ) eye " denotes the even part o/Todd .(x 1 ,x 2 ) under (x 1; x 2 ) i-» 
(-x l5 -x 2 ). 

Proof. From the definition of L A (S), we have 

L X (S) = L X ~\-S). 

Let A r ; := e 2ni ( ai ' r ) = Xo,i{r) for i = 1,2. Then we have 
2Todd* ven (x 1 ,x 2 ) =Todd £T (x 1 ,x 2 ) +Todd a (-x 1 ,-x 2 ) 

=sn(^ i ^)+f)+sfi(^ i (-^)-|) 

yer 6 i=l yer 6 i=l 

|r CT | 

=2L a {x 1 ,x 2 ) + —x 1 x 2 
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This finishes the proof. 

□ 

We say two cones cjj and cr 2 are similar if 

AA — — A~ 

for A e GL 2 (Z). In this case, A induces an isomorphism of M a in M a2 , 
which descends to isomorphism of T a in Y a2 . Since this isomorphism 
takes the lattice generators of a 1 to those of cr 2 , the two characters are 
preserved. A priori the Todd series of two cones coincide. 

Proposition 4.3. For two similar cones a and %, we have 

Todd cr (x l3 x 2 ) = Todd T (x l3 x 2 ). 

Proof. Clear. □ 

One should be aware that this is not the similarity of the matrices in 
linear algebra. cr(v l5 v 2 ) and cr(v 2 , v a ) are not similar in general. 

4.3. Dual cone and its lattice. LetiV := Hom(M,Z) be the dual lattice 
of M. N is a lattice in the vector space iV R := N (g) M. Using the standard 
inner product (x,y) we will often identify N and M. Associated to a 
lattice cone cr in M, its dual cone a is defined as 

d-:={yeiV R -0|(y,x) > 0} 

As the orientation is concerned, a is endowed with the orientation given 
by transpose of the matrix of a. Notice that a is again a lattice cone 
generated by two primitive lattice vectors inward and normal to cr. To 
d", there are two lattices naturally associated. M & is the sublattice of 
M = N generated by the primitive lattice vectors of cr. Note that this 
coincides with the definition of M a in 14.21 N d := Hom(M (J , Z) is a lattice 
in N R generated by dual vectors a x , a 2 to v 1 ,v 2 if cr = cr(v 1 ,v 2 ) (ie. 
(a u Vj\ = 5 i; ). These are related by the following inclusion relation: 

4.4. Euler-Maclaurin formula for 2-d cones. Let L k ' x (S) be the trun- 
cation of I A (S) below degree k: 

L k '\S) = (i + — )S + Q 2jA (0)S 2 + • • • + Q K MS k - 
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For a cone a, define a weight function as 

1 x g int(a), 
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wt K a sw (x] 



1/2 xeda-O, 
1/4 x = 0, 
otherwise 

Then the Euler-MacLaurin formula (with remainder) d la Karshon-Sternberg- 
Weitsman in | |22 || applied to / (x) evaluates the summation of the values 
of f(x) on the lattice points of cr w.r.t. wt™ '. 

Theorem 4.4 rKarshon-Sternberg-Weitsman||22||'). Let a be a lattice cone 
such that 

a = ct(u 1 ,u 2 ) for u ; e Z 2 . 

Let a ; be a dual basis of u fie. l^a^u^ = 5^) and X r>i := e 2jT! ( r ' a <) for 
i,j = 1,2. Let 

cr(fr) = {x gM 2 |(i,u ; ) ^ -fr ; , /or i = 1,2} 

/or h = (h 1 ,h 2 ) e M 2 . Suppose f is a smooth integrable function on cr(h) 
for some h e (K + ) 2 . Assume further that f is rapidly decaying toward the 
infinity of cr. Then we have 

2 wtf*(x)/(x)=2 f f(x)dx +Rl(f). 



xetjn 

where 

Kif) 



2 r°o rco 

|r*l r tF d Jo Jo r> 2 r> 
i r°° r°° 

lUl^r. Jo Jo 

i r°° r°° 

I r *l r trt Jo Jo 

In the above theorem, the summation is weighted with wt™. Now 
we come to the same theorem with the weight changed to wt a in (12. ID 
for our purpose. 

Corollary 4.5. With the same notations as in Thm \4.4\ we have that 
2 wt a (x)f(x) = \ 2 (UL k Md hi )-jd hl d h2 )] f f(x)dx +B°(f). 

xeanZ 2 ^ i=l 4 ^CO 
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Proof. This comes from the following observation: 

f(0,0) = d hi d h \ f{a l x 1 + a 2 x 2 )dx 1 d: 

J—h 2 J— hi 

= \r*\d hl d h2 f f(y)dy 

Jalh) h=0 



and 



2 wUx)f(x)= 2 wt K a sw (x)f(x)-lf(0,0). 



□ 



4.5. Asymptotic expansion. For A = e 2m i an iV-th root of 1 and k ^ 
2, from the periodicity one sees that Q k ^(x) is a bounded continuous 
function. Since we are assuming that / decays rapidly toward x,y — > 
+oo, we have 

J o J o ^ a ^ 2 /( a i x i + a 2 ^2)d^id^2 

is bounded for any i, j ^ 0. 

Thus the summation of the values of / (tx) over x running over lattice 
points of a is expressed in series occurring in Euler-MacLaurin formula 
involving the L -series. 

Theorem 4.6. Under the same assumption on a as in Thm. \4.4\ we have 
the asymptotic expansion for t e R + and t —*■ 

V wt a (x)f(tx) ~ (Lz(d h ,d h2 ) - ^-d hl d h2 ) o f f(tx)dx 
Proo/ We note that 

f(tx)dx = t~ 2 f(x)dx. 

a(h) Ja(th) 

Applying Thm. 14.41 for t e R + , we have 

2 wta (x)/(tx)=2 (t- 2 n^(t^)-A^) f /(*)d*| +jw)(t), 

xecrnZ 2 rer> i=l ! Ja(h) 1 
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with 



1 r°° r°° 

i r 00 r 00 

+ ^fS(- 1 ) fc_ J o J o Q fc) A r , 2 (x 2 )^I fe ^, 1 (_^ i )/(a 1 tx 1 + a 2 tx 2 )dx 1 d 

+ ^f2(- 1 ) fc_1 J J Qk* r Jxi)^K-d*Kf(aitXi + a2t*2)dx 1 d 
iorr = r d . 

Let us change the variables: y i = tx t for i = 1,2. Then from the 
boundedness of Qk,x ri > as /(^jj) decays rapidly, one can easily see for 
example a summand of R£ (/)(£) in the second line belongs to 0(t fc_1 ) 
at 0: 

J J Qk,^ 2 (x 2 )d k X2 L k ^^d Xi )f(a 1 tx 1 + a 2 tx 2 )dx 1 dx 2 

= t ^Jo Jo Q ^(T ) ^ lMrA(_t ^ )/(aiyi + a23r2)dyidy2 



= o(t fc - 1 ) 

Similarly one can show the rest belongs to 0(t fc_1 ). Thus we conclude 
that 

Rl(f)(t) = 0(t k - 1 ). 

□ 

The asymptotic expansion in theorem l4.6| can be rewritten using todd 
power series from lemma [4721 

Theorem 4.7. For t e 1R + and t — > 0, we have the asymptotic expansion 

2 wt ff (x)/(tx)~(Tbddn^A)-I^A) o f /( tx ) dx 

-a V Z / Jcr(h) 



fi=0 



4.6. Evaluation of zeta values. Now we apply this to the exponential 
series associated to the partial zeta function of an ideal. We saw that 
a partial zeta function can be identified with a zeta function of a cone 
cr w.r.t. a quadratic form Q. We are going to apply Zagier's theorem 
to evaluate the special values of £(b,s) at non positive integers. We 
apply Thm. 14.71 to obtain the asymptotic expansion of the exponential 
series. We take f(tx) = e -0-(t ll2 x 1 ,t 1 ' 2 x 2 ) = e -Q(x u x 2 )t ^ Then we b tam t ^ e 
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following asymptotic expansion: 



2 H* a e^~ 

ietxnM 

Theorem 4.8 f Garoufalidis-Pommersheim ||12|| 1 . For n > 0, we have 
C(b,-n) = 

(-l)»nl{lbddf +2 \c\,d h2 ) 5 n Ad hl d h \ o f e-^^)dx a dx 2 
where 

\ ifn = 0, 



(h,h 2 )=0 



5 n,0 



otherwise 



Remark 4.9. Jn ||12||. they used the exact Euler-Maclaurin formula of 
Brion-Vergne([4]) to a cone over a polytope. We apply the version with 
remainder term due to Karshon-Sternberg-Weitsman in [22 J. Since this is 
a formula with remainder, we have a direct method to evaluate the zeta 
values at nonpositive integers. 

5. Additivity of Todd series and cone decomposition 

In this section, we recall some technic used in | |12 || concerning addi- 
tive decomposition of the Todd series w.r.t. cone decomposition. 

Todd series does not allow decomposition in its original shape until we 
normalize. The normalized todd power series S a for a cone cr is defined 
as follows: 

^^-^^ de^!)^ 70 ^^^- 

One should note that different choice of the orientation of the same 
underlying cone yields the opposite sign in the normalized Todd series 
and interchanges the two variables. This is contrary to the original Todd 
series case, where the similarity class is determined by the sign. 

Let V; e 1R 2 for i = 1, 2, 3 be pairwise linearly independent primitive 
lattice vectors in a half-plane. An ordered pair (v i; v ; ) for i ^ j deter- 
mines a lattice cone cr^ = a^v^Vj) with orientation. 

In this case, we write formally 



Then we have the following: 
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Theorem 5.1 (Garoufalidis-Pommersheim [jl2[p. For i = 1, ••• ,r + 1, 

let Vj be pairwise linearly independent lattice points in a half plane of M 2 . 
We define cones 

cj; := o-i(v £ ,v i+1 ), o := cr(v l5 v r+1 ) 



cr = CJi + cr, + • • • + cr 



Thus 
Then 

r 

S a (x 1 ,x 2 ) =J]s ai (A-*A (7 (x 1 ,x 2 y). 

i=l 

In particular, if every a ( is nonsingular(i.e. det(A (T ) = +1) for i 
l,2,...,r, 

r 

S a (x l3 x 2 ) = 2det(A £7i )F(A a 1 A t7 (x 1 ,x 2 ) t ), 



i 



where F{x u x 2 ) } _ c t ^ „. 

Proof. See Thm. 2 in [[28]]. □ 

Remark 5.2. Abusing the notation, we denote cr(v 2 , v x ) by — cr(v l3 v 2 ). 
Actuary by definition of Todd power series of cone above, we easily find 
that 

Todd cr (x lJ x 2 ) = Todd_ tJ (x 2 ,x 1 ). 
The matrix A' 1 represents the linear transformation v 1 >-> e x , v 2 >— ► e 2 . 

So we have AZ^ = K At- Let A^ 1 = (^J f or two row vectors w i> 

,-i _ / W 2 



w 2 . Then A a = I I. Therefore, 

Todd a (A; 1 (x l3 x 2 ) f ) = Todd a ( (w 1; (x a , x 2 ) ) , (w 2 , (x a , x 2 ) ) ) 

(5.1) = Todd_ CT ((w 2 , (x l3 x 2 )) , (w 1; (x l3 x 2 ))) 

= Todd_ ff (A: 1 cr (x 1 ,x 2 ) t ). 

Thus one can see easily that for the additivity theorem to hold the orienta- 
tion of <j does not make any problem. 

6. Cone decomposition and Continued fraction 

In this section, we will decompose the cone a(b _1 ) into nonsingu- 
lar cones. This decomposition follows directly the decomposition of the 
fundamental cone in the totally positive quadrant of Minkowski space 
under the action of the totally positive unit group. This is fairly stan- 
dard fact related to desingularization of a cusp of the Hilbert modular 



20 



BYUNGHEUP JUN AND JUNGYUN LEE 



l 2i 


B = l(1) 

^ . * B_ 2 = a s _ l B_ 1 + B 




^\ „ / \ Bo = a.„_?B_? + Bi 
^B_ x = i(a>) 



Figure 2. B ; and the continued fraction [[a , a 1 ,..., a r _ 1 ]] 

surface of the real quadratic field considered. It is described in terms 
of the (minus) continued fraction expansion of the reduced basis of b" 1 
so that the desingularization of the lattice cone o"(b _1 ) in the sense of 
toric geometry follows(cf. ||11||. II14II ). We are going to apply Thm. 15.11 
to obtain explicit formula of the zeta values using the terms of the posi- 
tive continued fraction. One should note that our expression is differed 
from | |12 || in that we use the positive continued fraction instead of the 
negative one. 

In general, there are many other decompositions possible for a singu- 
lar cone. But those lattice cones arising from a singular cone but for a 
totally real field and the action of the totally positive units has a decom- 
position after the shape of its Klein polyhedron which is a geometric real- 
ization of a continued fraction. In 2 dimension, this appears as follows: 
In each quadrant of K R , we take the convex hull of b 1 = [1, co] and 
union the polygonal hulls. This is the Klein polyhedra of the ideal lattice 
b _1 . One can further assume co is a reduced basis :co > 1, — 1 < co' < 
for co e K. Then co has purely periodic positive continued fraction ex- 
pansion 

co = [[a ,a ir -- ,a r _i]]. 

Let {B 2i } (resp. {B 2i+1 }) be the vertices of the convex hull of t(b _1 ) 
in the 1st (resp. the 4th) quadrant of 1R 2 with B Q = t(l),B_ 1 = t(co) 
and x(B i ) < x(B i _ 1 ), where x(-) is taking the 1st coordinate. These B t 
arising as the vertices of the Klein polyhedron should not be confused 
with the Bernoulli numbers. 
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Let I be the even period of the continued fraction expansion of co(ie. 
I = r (resp. 2r) for even r (resp. for odd r). 

B ; satisfies a periodic recursive relation read from the continued frac- 
tion of co (cf. II14II ): 

(6.1) B i _ l = a i B i + B i+1 

Since a successive pair B h B i+1 is a basis of the lattice i(b _1 ) in K R , this 
yields a change of basis 




After successive change of basis, we have 

(6. 2 , ».,n)-i.^(v ;)(v J) 

Let otj, /3 f be the coordinate of w.r.t. the basis {B_ 1 ,B }: 

Note that the coulumn vector (a i ,j3 f ) t is equal to the 1st column of 

the matrix (Y j) (V j) " ■ ( V j) in (El. 
As B; are primitive, so is (a i3 /S ; ) in M = 1? . 

In the following, the totally positive fundamental lemma is identified: 
Lemma 6.1. Let e be the totally positive fundamental unit of K. Then 

t(e) = B_ t = a i _ l B_ l + ft_iB . 

Proof. Let e K > 1 be the fundamental unit of K. Then for the period r of 
continued fraction expansion of co, we have 

See p.40 of QT4D for detail. Since the totally positive unit e is either e K 
or e 2 K according to the sign of L 2 (e K ), we then obtain that t(e) = B_ e . □ 

Recall that we associated a lattice cone cr(b _1 ) in E 2 to b _1 in Sec. [21 

(6.3) cx(b- 1 ):=c7((0,l),(a,_ 1 ,^_ 1 )). 

This corresponds to the cone bounded by t(l) and i(e) in K m . For the 
rest of this section, only the cone cj(b _1 ) is need to consider to compute 
the zeta values. So we will write simply cr instead of cr(b _1 ). 

Lemma 6.2. Let a = cr((0, 1), (a, jS)) be a lattice cone where a,/3 are 
relatively prime positive integers. Then a the dual cone of a is similar to 

T = T((0,-l),(a,/3)). 
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Proof. It is easy to see the dual cone a has primitive basis ((1, 0), (— j8, a)). 
See Fig. [31 Since the rotation by —90 degree belongs to SL 2 {Z), we have 
the desired similarity of the cones. □ 

After Prop. 14.31 and Lemma l6T2l for b as before, we have 

o- ~ T := T((0,-l),(a f _ 1 ,^_ 1 ) 

thus 

Let v_j = (0, 1), v = (1, 0) and for 1 < i < I - I, 

for a ; ,/3j defined as in Eq. (l6.2D . Notice that v t corresponds to B i+1 and 
v_ l3 v are the two standard basis of M. Then the decomposition of a 
yields that of cr: 

Proposition 6.3. With above notations, let cr' := cr'(— v_ l3 v ) and 

(j. ■= (JiiVi^Vi), 

for i ^ 0. Then we have 

a ~ t := t(v_x, v £ ) = cr' + (Tj + tr 2 + cj 3 + • • • + a^. 
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Thus we have 

i-i 

S & ( Xl ,x 2 ) = F(A-}A T ( Xl ,x 2 Y) + Y i {-^ i F{KlMxi,x 2 ) t ). 

i=l 

Proof. First, one should notice that 

T = a' Q + p 

where p = p(v , v £ _ 1 )(See Fig. [3]). Note also 

det(/g = det (jjj-j a A = fta,_ a aA^ = 

Hence the decomposition of p into nonsingular cones a i 

P =cr 1 + cr 2 + --- + cr £ _ 1 

finishes the proof. 

Lemma 6.4. For -1 < i < I - 1, let 

Mi := (-l) i+1 ((A^-i - OiPi-i)x 2 + a i x i)- 

Then we have 

/ 1-2 1 

Todd*(x l3 x 2 ) = a t _ lXl x 2 J] (-l) l F(M i} M i+1 ) + — - 

\i=-i 1 e 

where F(x 1 ,x 2 ) = . 

Proof. After simple computation, we obtain that 

A-^A T (x l3 x 2 ) t = (M i+l3 M i ) t . 

Since F(-x l3 x 2 ) +F(x l3 x 2 ) = 73^, we have 

F(C^ 1 T(x lJ x 2 ) t )+F(C- 1 T(x l3 x 2 ) t ) 
=F(x a - &_iX 23 cL t _ x x 2 ) + F(-(x a - &_iX 2 ), a £ _jx 2 ) 
1 



1 — e- a i-i x 2 

Note 



det (-°i £;)-«<- 

If we apply the above to Prop. I6.3L we complete the proof. 

Let Todd (T (x 1 , x 2 )^ be the degree n homogeneous part of Todd^f 



&2k+2 

(2k + 2)\ 
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Proposition 6.5. Let 

2fc + l R R 

L k (X,Y)= V -i - 2fe+2 ~ i x^Y 2 ^ 1 , 

fcV ; £J i! (2fc + 2-i)! 
R fc (x,y) =X 2fc +X 2fc - 1 7 + --- + 7 2fc . 

Then we have 

Todd ( ,(x 1 ,x 2 )( 2fc+2 ) = 

(£-2 l-l 
2 (- 1 ) i - L fe( M i+i 3 M >i^2 + 2(-l) i ^-i^(M i _ 2 ,M ; )x 1 x 2 
i=-l i=l 

+ (2 ^ + 2), (- y i< +1 + ^MS 1 ) + S^a^x^. 
Proo/ From Prop. |6.4[ we find that 

Todd^x,)^ 2 ) = a e _ lXl x 2 2 (-^^(M^M^O^+X! ^1/^ . 

;=-i 1 e 

We have 

F(M i; M i+1 )( 2fc ) 

2fc+l R p R /m 2Ic+1 M 2k+1 
= y ^m_f2fc+2_m_ Mm _i M 2fc- m +l ^fc+2 / Mj + l 

^ m! (2k + 2-m)! i+1 ; (2k + 2)! y M i+1 M f 

=L fc M i+1 ,M ; + - , +2 — ! + — 

fcV !+1 i7 (2fc + 2)! I M i+1 M f 

and 

^-l X 2 (2fc+1) 2fc + l,.2fc + l , 1 „ 

U l~l 2 ~ U k,Q o u -(.-l x 2> 



l-e-«s-i*2 (2fc + l)! ^ 2 *' u 2 

as B 2fc+i = for fc> 0. 

Moreover we also have the following: 

1-2 

2 (-iy(M-\M 2fc+1 + Mf^Mr 1 ) 
(6.4) i=_1 

M 2fc+1 Mf+ J .M 2 * +1 -M 2fc+1 

As 

'(o_ 1) j8_ 1 ) = (0,l) 

a i+l = a «-i-l a L + a i-l 

A+i = <k-i-iPt + Pi-i 
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we have 

M_ a = ^_ia £ _ x x 2 - a_ a (-x 2 + /3 £ _!X 2 ) = a^Xa, 
^-i = &-i^-i* 2 - + ^-1*2) = a<-i*i- 

and 

Therefore Q6.4D is equal to 

M n 2fc+1 Mf+ X ^ M 2fc 2 +1 -M 2fc+1 

- — — + — — + Yx-iya^— ■ 

(6 ' 5) M 2k+1 M 2k+1 l - 1 

' ° + T ^- + 2(-l)V 1 ^(M i _ 2 ,M i ) 



Thus we finally complete proof. □ 

7. Special values of zeta function 

Now we are going to evaluate the values of £(s, b) at non-positive 
integers using the expression of the degree n homogeneous part of the 
Todd series made in the previous section. We suppose b is an integral 
ideal normalized as in the previous section so that b 1 = [1, co] for 

co = [[a ,a l3 --- ,a r _ 1 ]]. 

and e > 1 denotes the totally positive fundamental unit of K. Let I be 
the even period of continued fraction expansion of co. 
(a u /3 ; ) for i = 1, 2, ■ ■ ■ , I — 1, and 

(a_ 2 ,jS_ 2 ) = (l,-a ), (a_ a ,j8_ a ) = (0,1), (a ^o) = (1,0) 

are primitive lattice vectors in M. Note that (a ; , /3J corresponds to B_ i+1 
inK R . 

Q(x 1 ,x 2 ) := Nityfaa) + x 2 )(x 1 co' + x 2 ). 
Then the partial zeta function £(s, b) is expressed as (See Prop J2.in 

^ = 2, Q(Z)' • 

where 

C7(b- 1 ) = a((0,l),(a,_ 1 ,ft_ 1 )). 

In Thm j4.8L the partial zeta value is written using the Todd differen- 
tial operator of the cone a dual to cr = cr(b _1 ). We apply the additivity 
of the Todd seriesfProp J6.50 after the cone decomposition of a occur- 
ring in the continued fraction of co to this expression. Then we obtain 
the following expression of the partial zeta value: 
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(7.1) Z(-k,b) = (-l) k k\(^ + &)o ( e- Q( - Xl ' X2) d Xl dx 

JaQi) 



h=Q 



where 



(7.2) 



and 
(7.3) 



<£ : = J] (~l)^ fc (M i+1 ,M ; )(^ i ,4 2 )^-i4 1 ^ 



t=-i 

l-i 



- ^J^L S(-l) i a<-^ fc (M f _ 2> M i )(4 i> g h2 )a < _ 1 g hi g h 



(2k + 2)!^ 



7^^ (-a f R fe (M_ 2j M )(4 lJ 4j^-i4 1 4 2 -4 1 M 2fc+1 (4 1 ,4 2 ) 

+ 4 2 mS 1 (4 1 -^ 2 ))' 

In (I7.1D , as the differential operators are linear, this expression can be 
evaluated one by one. Later in Secj9]we will prove that the part of (I7.1D 
involving 52 vanishes: 

(7.4) & o I e - Q(xi ' X2) dx 1 dx 2 = 

JaQi) h =° 

Then it remains only to evaluate the part involving i£ . First, we need 
to rewrite the integral in another coordinate (y l3 y 2 ) such that (x 1 ,x 2 ) = 
(ct^_ 1 y 2 ,^_ 1 y 2 + yO" So we have cr(h) in the new coordinate: 
a(h)=a(h l ,h 2 ) 

= {yi^i + y 2 v 2 |(y 1 v 1 +y 2 v 2 ,u 1 ) ^ -h 1 ,{y 1 v 1 + y 2 v 2 ,u 2 ) ^ -h 2 } 

h x h 2 

= {{at-iy2,Pt-iy2 + yi)\xi > ,x 2 > }. 

a i-\ a i-\ 

In the new coordinate (y l5 y 2 ) the integral becomes 
(7.5) 

c r 00 r 00 

e -Q.{x lJ c z ) d d = e -Q^t-iy2,Pt-iy2+yi)d y dy 

JaQi) ' " J-JL- J__*L. 

v ' a i-i a t-i 



This integral applied by a^d^d^ is 
(7.6) a f -AA f e- Q ^^)dx!dx 2 = e 

JaQi) 



-N(b)N(-^-e+-^-) 
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h=0 



The above simplifies (I7.1D quite much assuming the vanishing of (I7.4D : 
(7.7) 

1-2 

((-Kb) = (-l)W J] L k (M i+l ,Mi)(d hi ,d h2 ) 
t=-i 

+ (^1)! 2 1- 1 )VA(M,,,M i )(4„y) oe- , " h,, "^T' t ^ , 
Lemma 7.1. Let A; = a t co + /3 ; . For — 1 ^ m, Z < I — 1, we have 

M i (4 1 ,4J ! Mj4 1 ,4J^- W(b)W( ^ e+ ^ ) |^ = 

/or M^x^Xa) = (-l) i+1 ((A«<-i - «iA-i)^2 + 

Proo/ For simplicity, let c ; = (— l) I+1 (i3 t ' a £-i — a iPe-i) andd ; = (— l) !+1 ctj. 
By internal change of coordinate (h 1 ,h 2 ) >-> (ah 1 + cfr 2 , fe/?! + dh 2 ), we 
have 

(a4 1 + fc4J i (c4 1 +d4J j /(^i> 2 )l^o = ^ 1 ^ 2 /(ah 1 +ch 2 ,bh 1 +dh 2 )|^ . 
Thus 

-N(b)N( h2 e+ hl ) 

HdicX+cAJid^ + c^Xyoe U + ^ J | h=0 
We note that 

p l . 1 d i + c i = {-l) i+1 p i a l _ 1 . 

Since e = a t _ 1 co + /3^_ 1; we have 

(d i h 1 + d m h 2 )e + c,h 1 + c m h 2 

= ( d iPi-i + c i) h i + ( d mPi-i + c m )h 2 + (dihx + d m h 2 )a £ _ 1 co 

=a t _ 1 (((-lY+'aA + (-l) m+1 a m h 2 )co + (-l) 1 * 1 ^ + (-l) m+1 AA 

=a,_ 1 ((-l)'+%h 1 + (-l) m+1 A m fr 2 ) 

□ 

We note that 

iV(b)N((-l)' +1 Ah 1 + (-l) m+1 A m h 2 ) 

= Q((-l) i+1 a i h 1 + (-l) m+1 a m fr 2 , + (-l) m+1 ^ 2 ), 
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for a binary quadratic form Q(x,y) with degree 2 and i, j with i +j = 2k, 
we have 



hi n 2 



h=0 k\ 



\fLdh\dh } 2 Q.{h lt h 2 ) 



h=0 



Thus, from (17.61) (I7.7D and Lemma 17.11 we have finished the proof of 
our second main theorem (Thm. 11.21) . 

i-i 



! = 

t-1 



h=0 



+ 



Blk+2 

(2fcTY)!^ 



£(-1) V^nKAjQl^Wh + a^fi^h, + fth 2 ) 



h=0 



Remark 7.2. We have obtained a polynomial expression of the zeta value 
in variables a^fii and the coefficients of the quadratic form Q(x,y). For 
the polynomial expression, it is important to show the vanishing Q7.4D . In 
Sec. [9] there appear a , ct l _ 1 in the denominator of the vanishing expres- 
sion involving the 8%-operator. This is a crucial ingredient of the Kummer 
congruence and the corresponding p-adic zeta function. 

8. Computation of £(— k, b) for k = 0,1 and 2 

In this section, we evaluate the zeta values £(— k, b) explicitly for small 
n. We express the values in terms of the continued fraction expansion 
[[a , a l3 . . . , a^_i]] of the reduced basis co of b _1 . 

8.1. k=0. For k = 0, the zeta value is already known by C. Meyer(fl27l]) 
in terms of negative continued fraction. Using the plus-to-minus conver- 
sion formula of continued fraction 

(8.1) 

5 = co + 1 = [[a , a 1 ,..., a^]] + 1 
= ((a + 2,2,...,2,a 2 + 2,2,...,2,a 4 + 2,...,a^ 2 + 2,2,...,2)) 
= (Oo> b n---> fe m -i)) 

one obtains the result in positive continued fraction. 

In our approach, we begin with the expression using positive contin- 
ued fraction as a special case of Thm J1.2t 

«o,b) = ^2(-i)V 1 

z i=o 
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Since B 2 = 1/6 and I is the even period of the continued fraction, this 
reduces to 



1 l ~ l 

(8.2) a0,b) = -2(-l) 1 « 1 

iZ ;=o 

Via (I8.1D . one recovers the result of Meyer: 

2 m— 1 

c(o,b) = 3), 

iZ ;=o 

where b ; is the i-th term of the negative continued fraction. 

Remark 8.1. Note that using the positive continued fraction, we have an 
alternating sum for the zeta value. Consequently, one sees directly the van- 
ishing of £(0, b) when the actual period of the positive continued fraction 
of co is odd (equivalently, if the fundamental unit is not totally positive). 

8.2. k=land2. For Q(x l5 x 2 ) = N(b)N(x 1 co + x 2 ), let L h M t and N t be 
defined as in 

Q(a i h 1 - a(_ a h 2 , PA - &_ x h 2 ) = + M f M2 + tyh 2 . 
Similarly, L i3 M, and JV f are defined as follows: 

Q(a i _ 2 h 1 + a ; h 2 , ^_ 2 h a + /3 ; h 2 ) = Ljh^ + + iV ; fr 2 

Then the special value at s = 1 is computed out as follows: 

b) = ^(-l)'- 1 ^ + 5l 2(-l) ! a,_ l (2L ! +M l + 2N t ) 

i=0 ^ 4 " i=0 

1 £_1 

= ^ ( 5M > + a ^-i( 2l t + + 2 ^)) 

/AV i=0 

Similarly for s = —2, 

a-2,b) = 
i 

-----^ ^(-l) 1 (21M ; (iV; + L t ) + 2a t _, (6Lf + + M 2 + 2L ; iV; + + 6JV 2 )) 

i=0 

This should be compared with the expression obtained using nega- 
tive continued fraction in | |12 || and also H34II . They considered the zeta 
function of the following quadratic form in view of negative continued 
fraction: 

Q'(x l3 x 2 ):=JV(b)JV(x 1 5 + x 2 ) 

for 5 = co + l(See (18. ID for negative continued fraction). Let A t be the 
lattice points of the component of the Klein polyhedron of b 1 in the 1st 
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quadrant with normalization: A = 1, A_ 1 = 5 and the 1st coordinate of 
A; increasing according to i. Then we associate a lattice vector (p k ,q k ) 
to A k for A k = — Pf c A_ 1 + <2fcA)- p k and q k are obtain from the reduced 
fraction of the truncation after k the of the negative continued fraction 
5 = ((b , b 1 ,...,b m )): 

— = (b ,... 3 b fc _i) 
Pk 

(This is the last line of pp.18 of [12], where — should be corrected to — 

ik Pk 

as we just wrote above). Similarly, L' M ; ', N! and L i ,M t ,N t are defined 
as the coefficients of quadratic forms: 

Q'i-Pi-ih Pih 2 , qi _ r h x + qi h 2 ) = L'h\ + M[h x h 2 + N[h\ 

and 

Q'i-Pi-iK p i+1 h 2 , q^h, + q i+1 h 2 ) = L[h\ + M/h^ + N!h 2 2 . 
In this setting, Garouf alidis-Pomersheim ( | 1 1 2|| 1 obtained: 

m-1 

720 

and 



cc-i, b) = ^ 2 ( 5M / + b i(~ 2L i + - 2i ^)) 

i=0 



2 m— 1 

a " 2 ' b)= iii2o2(- 2iM ! '( i ;+ iv ! / ) 

i=0 

+ 2b ; (6L' 2 - 3L'M' + 2L'.N'. + M' 2 - 3M'.N'. + 6N' 2 )). 
i ^ i ii ii ii i it 

It should be also compared with Zagier's result(eg. for k = 1) in [34]: 

m-1 

720 



b) = ^ 2 (-2N i bf + 3M ; bf - 6L ; b £ + 5M()) 

i=0 



Remark 8.2. if we use the formula for zeta values using negative contin- 
ued fractions as is made by Garoufalidis-Pommersheim and Zagier, one can 
still obtain polynomial behavior in a family similar to Sec. \T0\after the uni- 
formity of the negative continued fractions in the family. But as is known, 
there is hardly a direct arithmetic property(eg. regulator) associated to 
the negative continued fractions. Actually in ||20||. ||21||. the formula for 
negative continued fractions, which is developed by Yamamoto( [33]) and 
Zagierf []3"4l , ||35||). is used after conversion of positive continued fraction 
into negative one. The formula using positive continued fraction simplifies 
this unnecessary step and justifies the reason of our earlier results. 



special values of zeta functions 

9. Vanishing Part 
Now, it remains to show the vanishing of (I7.4D 



31 



St q I e- Q(xl - X2) dx 1 dx 2 



0. 



h=0 



This part is crucial in expressing the zeta values at nonpositive in- 
tegers as polynomials of its argument coming from terms of continued 
fractions and the coefficients of quadratic forms. The vanishing has been 
already observed in related works by Zagier([34]) and Garoufalidis- 
PommersheimC [fl"2~H ) in different settings. In [12], only the vanishing 
is mentioned without clear proof. In this section, we will recycle some 
notions and ideas from II34II . 

It suffices to show the vanishing of the following, which equals the 
above up to multiplication by a constant. 



(9.1) 



-a i R k (M_ 2 ,M Q )(d hl ,d h2 )a l _ 1 8 hl 8 h2 - d h M^ +1 (d hi ,d h . 



+d h2 MfZt 1 (8 hi ,d h2 )) o f e-o^dx^x 



h=0 



As M = I3(_iX 2 — x 1 and M t _ 2 = x 2 — a e _ 2 x 1 , we have 



(9.2) 



2xf + 2 + 



\ i 2fc+2-i 



Applying d^ +2 to ( 17751 ), we obtain 



(9.3) 



a t-i 

i f°° 

a t-l Jo 



,-N(b)N(ey 2 + yi ) 



a l-i a l-i 



jV(b)JV(e^2 !a-) 



dyidy 2 
dy 2 



h=0 



h 1= 
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IfwewriteP(x 1 ,x 2 ) = ^P-(xl + (e + e')x 1 x 2 + x^), using (l9T2l)-(l9T3lh 

a t-i 

one can simplify the 2nd half of (I9.1D : 
(9.4) 

- d h M (c\, d h2 ) 2k+l + d h M e _ 2 (c\, d h2 ) 2k+1 ) o f e-^^dx.dx 



+ 



a «-l Jo 



02k+l e -P(- Xl ,x 2 ) 



dx 



From M„ 2 = ( a o a «-i + fle-i) x 2 ~ x i an d = Pt-i x 2 ~ x i> we have 
(9.5) 



2k+l 



R k (M_ 2 ,M )= ^(-l)' 



+i 



i=0 



2lc + l\ (VH + ft-il'-ft'-i J _ 1 2k+1 _ 



X 2 X l 



From (I9.4D and (19.5D , we obtain the following lemma: 
Lemma 9.1. Let P(x l3 x 2 ) = -^-(x^ + (e + e')x l x 2 + x^). Then we have 

- a^ fc (M_ 2 ,M )(4 i5 4 2 )a^_ a 4 1 4 2 - d h M* k+1 (d hi , d hz ) 

+ d h Mf n - + 2 \d hi ,d h2 )) o f e-Q^dx.dx 

dx 2 



(ao«£-i + A 



fi=0 



Jo 



fl2k+l e -P(- Xl ,x 2 ) 
x i 



Xi=0 



Lemma 9.2. For the totally positive fundamental unit e > 1, we have 

e + e' = a Q a l _ 1 + p t _ x + a e _ 2 . 

Proof. We note that 



5 := = [[ai-i,CLi-2,-'- , a o]]- 

CO 



Thus 



a r 



1 
And we have 



a l-2 1 

1 



a 1 \ f 5 \ _ a^ + q^i 
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Finally we have 



co + CO 



, a i ~ Pi-i 



a, 



-l 



Thus 



e + e' = a^co + co') + 2^_ a = a g + p t _ y 
e + e' - /3 t _ 1 - a e _ 2 



at 



□ 



From Lemma |9T2l if we let a Q a t _ l + = —a, a e _ 2 = — b and = 

a t-i 

A then we find that e + e' = — (a + b). Hence one can rewrite Lemma [97T| 

as follows: 

(9.6) 

2k 



i=0 ^ ' 



fi=0 



+ 2 | ^k+i e -p(- Xl ,x 2 ) 



dx 



Xj^O 



i=0 V / 



h=0 



+ 2 ( 52fc + l e -A(x2 + (a+b)x 1 x 2 +xJ) 



dx 



x,=0 



Hence it remains to show vanishing of the right hand side of Q9.6D : 

2k 



h=0 



2 ( 2 i+i) ( ai+1 + bi+1 )^?"° rA( ^ (fl+i),,, '' 2+ft ' ) 

(*) i=0 V 1 / 

poo 

+ 2 ^2fc+l e -A(x2 + (a+b)x 1 x 2 +x2) dx 

Jo Xl -i=o 2 ' 

For the proof, we introduce f k (a,fi,y) and d rk (a, as follows: 

(2fc + l)!/ fc (o,j8,r) 



(9.7) 



(9.8) 



^2fc+l e -(ax^+^x a x 2 + r xf) 



n=o - 2y k+1 



2/c 



2\k 
21 



i=0 



These numbers are originally appeared in H34II . One should see that 
f k (a,p,y) is odd function w.r.t. /3: 

f k (a,l3,r)+f k (a,-p,y) = 0. 
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One can identify d^ 2 fc_i(ct 3 — j8,y) in the following expression: 

<9 ; ^ 2fc -ig-( ax i-^ x i x 2+r^2) 

(*l,x 2 )=(0,0) 

= ^p<^- £ (axj - Px lX2 + rxlf = ^pi!(2k iVA^a, -/3, r ). 
From this, one can rewrite the 1st line of (j*j) as 

2k 



2 (^i + i) ( a ' +1 + b i+1 )5^ h 2 i fc - t 'oe-^-( a+b ^ + ' l2 1 ) 



h=0 



(9-9) 1=0 

(-l) fc |V +1 + b i+1 
= y - T ^ L (2k + 1)! J] 7— j d U2k _i(A, -A(a + b),A) 

i=o + 

The 2nd line of © is, from the definition of f k (a, /3,y), 
(9.10) 

2 r^ + i e -A(xl + ( a+& )x 1 x 2+ x ! ) d = _ (2 fc + l)! A(AjA( " + b)jA) . 
Jo Xl -1=0 V ' A k+1 

Now, we are going to use an identity relating f k ( a, /3, y) and d; j2 fc-i( a ) — y) 
due to Zagier: 

Lemma 9.3 (Zagier(Prop. 4 of [34])). For a real number X, we have 
f k (a, 13, Y ) + f k ( Y , 2Xy 0, X 2 y A/3 + a) 

(—l) k 2k X i+1 

i=0 

If we put a = A, /3 = A(a + b), y = A and X = a (resp. X = b) into the 
above, we obtain 

f k (A,A(a + b),A)+f k (A,A{a b),A(-ab + 1)) 
(-1)" 

= 2 y —^A k+1 J]d i!2k _ i (A,-A(a + b),A) — 
i=o 1 

and 

/ fc (A,A(a + b),A)+f k (A,A(b a),A(-ab + 1)) 

= 2^ T rA t+1 ^d ij2t _ i (A,-A(a + b),A)r—r. 

i=o 1 + 1 

As / fc is odd function of its 2nd argument, summing the above two 
equations, we have 

f k (A,A(a + b),A) (-l) fc £. , A , tU , A y +1 + b i+ * 
^71 = — L^k-iiA-Aia + b),A) . 

i=0 
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This identifies (l9?9l) and (f97TOT) up to sign. 
Therefore we concludes the vanishing of fig) . 

10. Appllication: Polynomial behavior of zeta values at nonpositive 

integers in family 

Until now, we developed a way to compute the partial zeta values 
at nonpositive integers for a real quadratic field with a fixed ideal b 
via the shape of the continued fraction of a> for b" 1 = [1, co]. We will 
apply this method to certain families of real quadratic fields to prove 
the main theorem of this paper(Thm. II. ID . We deal with the same 
family of real quadratic fields with ideals fixed as in our earlier works 
(" II20H , ||21||"). In our previous works, the partial Hecke's L -values and the 
partial zeta values of a ray class ideal at s = are investigated for family 
of real quadratic fields. We showed that the values in the family is given 
by a quasi-polynomial in variable n which is the index of the family of 
real quadratic fields considered. If the conductor is trivial, so that we 
consider ideal classes, the values behave actually in a polynomial. This 
method was originally observed by Biro and has been main ingredient 
to solve class number problems of the real quadratic fields in the family 
without relying on the Riemann hypothesis (cf. QU, QSJ, [0], [ITU). 
Here we deal with the case when the conductor is trivial. Thus we have 
strict polynomial instead of quasi-polynomials. 

We generalize the result on the partial zeta values at s = to every 
nonpositive integer s when the conductor is trivial. This means we con- 
sider partial zeta function of ideal classes instead of ray classes. So the 
scope of partial zeta functions we consider here is narrower than the 
previous. But the same method must be applicable to ray class partial 
zeta functions. In this case by the same reason quasi-polynomials are 
appearing instead of polynomials to give the zeta values at a given non- 
positive integer for the same family of ideals. Again this will answer the 
same for the partial Hecke L-values at arbitrary non-positive integers. 

Recall the conditions on the family (K n , b n ) indexed by n e N for a 
subset N of N. b^ 1 = [1, co(n)] for a reduced element co(n) e K n and 

a>(n) = [[ a oO)> a iO)>-" > a r-i(n)]] 

for polynomials a t (x) e Z[x] and the quadratic form N(b n )(xco(n) + 
y)(xo)(n)' + y) associated with b n is expressed as 

b 1 (n)x 2 + b 2 {n)xy + b 3 (n)y 2 

for polynomials b f (x) e Z[x]. 
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Proof of Thm. \l.l\ Applying Thm ll.2l to the family considered, we have 

2(-l) i " 1 ^(5 fll ,^ 2 )Q(a i (n)h 1 - a^n)^, ft(n)h a - &_ a (n)h 2 ) fc 

! = 

(2fc + 2)!^ 

Since Q(— ) is a quadratic form and a ; (n), b^n), a i (n),^ i (n) are poly- 
nomials, it is clear that £(— fc, b„) is a polynomial in n. 
Notice that 

dega; > degft 

and 

dega; ^ dega^. 

Thus, the highest degree term comes from the summand with i = I — 1. 
Putting altogether, we obtain the denominator C k as well as the degree 
m = kC + D for the explicitly given C, D. □ 

Remark 10.1. One should notice the independence of n of the denominator 
C k of £(— fc, b n ). A priori this is invariant in the family. It is important to 
control the denominator to interpolate the associated p-adic zeta function 
from the values at negative integers(cf. Q9D, II10II and II24II ). 

For the rest of the paper, for a number field K, let us denote the ring 
of integers by K . 



Example 10.2. Consider the family (K n = Q(Vn 2 + 2), b n = K J. Then 

b- 1 = Kn = [1, co n ] 

for co n = [[2n,n]]. 
Then we have 



12 

r(-i b ) = -— + - 



ao,K) 

_ 2159n , 
25200 

68n 797n 3 , 689n 5 , 134n 7 2878n 9 



360 40 

2n rf_ 23n 5 

45 945 1890 

2159n . 137n 3 59^ . 3n' 



X(_o U \ _ 2159n , 13/n 3 _ 59n 5 , 3rf_ 
^ ' 11 ' 25200 25200 840 56 



T '( — A h \ — — /y/n _i_ bBVn I 



231 6930 1155 1155 10395 

11947883n , 29660563n 3 26073083n 5 7603n 7 145933n 9 . 351719H 11 



y(_n (_ \ _. 1194/883n . ZVbbUbban" ZbU/iUU6n- 7bU3n' 14bV33n" 

' "J " 7567560 22702680 5675670 2310 135135 135135 
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Example 

Kn . Then 

C(-3,b n ) 



10.3. Let K n = Q( Vl6n 4 + 32n 3 + 24n 2 + 12n + 3) and b n = 
b n 1 = °K n = [1, "J M co n = [[8n 2 + 8n + 2, 2n + 1]]. 

n , 2n 2 



_ J_ 

~~ 12 



C(-4,b n ) 



, n , 2n^ 

2 3 
_7_ _ 13n _ lln 2 
72 20 



45 



34n 4 
15 



104n 5 , 32n 



45 



/ ^ y to id to 45 

_ 503 . 2773n , 8473n 2 , 13009n 3 _ 6898n 4 _ 360n 5 _ 6208n 6 
2520 1260 945 945 945 7 105 

_ 3328n 7 , 25472n 8 . 18944n 9 . 4096n : 
315 945 945 945 



823 
840 



32553n 4 , 143188n 5 , 
525 1575 

147072n 10 , 4862464n n 
35 1575 



2707724n 6 

1575 
. 830464n 12 
525 



105 

! 10 

yto yto 945 

7762n _ 193469n 2 _ 309377n 3 _ 232553n 4 
525 2100 1050 " " 

. 5759672n 7 , 7377392n 8 , 7421248n 9 

1575 1575 1575 

, 249856n 13 , 32768n 14 
525 525 

106613 262407k . 1957759n 2 23147174n 3 203979376n 4 365417032n 5 
_ 11880 1540 1386 3465 10395 10395 

, 257724232n 6 _ 764543312n 7 _ 1238665888n 8 _ 3134586496n 9 _ 3314036224n 10 

10395 10395 3465 3465 2079 

_ 4177427456n n _ 880111616n 12 _ 2179907584n 13 _ 27426816n 14 
2079 495 2079 77 

_|_ 12058624n 17 



_ 86638592n 15 , 22151168n 16 
3465 



I 16777216n 18 

3465 ^ 693 ^ 945 10395 
_ 4617527 _ 52647823n _ 71296254191n 2 _ 22288517357n 3 _ 2248765926611n 4 

34398 17199 2270268 ' 115830 " 2837835 

_ 202240251208n 5 _ 1639280941052n 6 _ 7365379306328n 7 _ 45593045200n 8 

85995 315315 945945 51597 

, 124808351658752n 9 . 501230433622144n 10 . 1216530615292672n n . 421424974443008n 

2837835 2837835 2837835 567567 

. 2733316068964352n 13 , 301312337145856n 14 , 694067022135296n 15 , 416253219504128n 16 

2837835 315315 945945 945945 

, 84428067700736n 17 . 223080016510976n 18 . 13420435865600n 19 . 15443400065024n 20 

405405 2837835 567567 2837835 

, 481111834624n 21 , 185488900096n 22 
567567 2837835 
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